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Abstract 

A representation without explicit use of the isospin formalism is developed for the precise 
study of few-nucleon systems, and the advantages of the proposed approach are demonstrated. 
Using the example of three-nucleon systems with central exchange NN-interaction potentials of 
the general type, the complete equivalence is shown for the both approaches with and without 
the isospin formalism. The new systems of equations contain a less number of components as 
compared to the commonly used ones and are more suitable for precise study of few-nucleon 
systems, in particular, within variational approaches. Optimal variational schemes are developed 
with Gaussian basis, and the binding energies, r.m.s. radii, density distributions, formfactors, 
and pair correlation functions are calculated with a high and controlled precision for 3 H and 
3 He nuclei. 
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1 Introduction 

Noticeable progress in precise calculations of the few-nucleon binding energies within variational 
approaches , @ , @ stimulates the further development of optimal variational schemes and their 
possible refinement in order to achieve higher precision in calculations || of the main structural 
characteristics of the bound states of few-nucleon nuclei. The progress in this field enables one to 
obtain "exact" results for various structural characteristics of the systems and to put forward the 
problem of precise fitting of the nucieon-nucleon potentials universal for all the light nuclei, as well 
as to carry out accurate variational calculations for heavier nuclei. 

The achievement of a high accuracy in calculations of a few-particle quantum system is restricted 
both by the number of particles and the number of the wave function spatial components. At the 
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first stage of studies of a simple qualitative type, rather rough assumptions were used about nuclear 
interaction potentials and the wave function spin structure (see ||, 0> H), thus an account of all 
the wave function spatial components was considered to be not necessary. But precise calculations 
with controlled accuracy need to take into account all the factors contributing to the result of a 
given accuracy. Such a statement of the problem is faced with rather nontrivial difficulties even for 
comparatively simple three- and four-nucleon systems. In the standard commonly used scheme with 
the use of the isospin formalism, an additional symmetry of the wave function and the mixture of 
states with different total isospins due to a consistent account of the Coulomb interaction lead to 
the essential growth of the number of spatial components of the wave function and thus to a serious 
complication of calculations. As a result, to find the total wave function of the 3 He nucleus in the 
framework of the isospin formalism, one has to solve a system of six spatial equations, while in the 
case of 4 He nucleus one has already to deal with twelve equations. 

Here, we present our approach proposed in || without explicit use of the isospin representation 
and, using the example of three-nucleon systems, show the complete equivalence of such a more 
convenient approach to the traditional isospin formalism. It is shown for the wave functions in both 
approaches to be connected by certain relations, while all the physical observables calculated within 
both approaches are shown to coincide. The obtained much more simple systems of equations for the 
spatial components of the wave finctions of three nucleons in the bound state or scattering process 
(in the doublet or quartet spin states) enabled us to carry out precise calculations for three-nucleon 
nuclei and to analyze their structure. 



2 Equivalence of the representations with and without 
isospin 

In order to show the complete equivalence of the representation without explicit use of isospin 
(protons are not identical to neutrons) to the traditional isospin formalism, we recall the explicit 
form of equations for the three-nucleon system (2p, n) within the isospin formalism with common 
assumptions about the central exchange nuclear two-particle potential and the Coulomb interaction 
between protons. In bound state ( 3 He nucleus), the system (2p,n) has spin S = 1/2 and the isospin 
projection T 3 = 1/2, being in a mixed isospin state (T = 1/2 and T = 3/2). Thus, the total 
antisymmetric wave function can be represented in terms of standard spin-isospin functions and 
consists of six spatial components: 

* a = re + (v>r - + re + w - v'Ox^ . (i) 

After projecting the Schrodinger equation onto the basic spin-isospin functions, one has the system 
of six equations for the spatial components of the wave function: 
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where K is the kinetic energy of three nucleons and 



U'33) V>' ~ (Ui 3 ~ E/33) if"} 



(2) 



U' = V(r 12 ) + V(r 13 ) + V(r 23 ), U' = ^-(-V(r 13 ) + V(r 23 )), U" = -V(r 12 ) + ^(V(r 13 ) + V(r 23 )) 

(3) 

are the combinations of spatial components of the potentials (both nuclear and Coulomb ones) with 
a certain symmetry with respect to permutations of nucleons. The last two equations in (|2|) are 
present due to the mixing of different isospin states by Coulomb interaction. If we set U c = in (0), 
we get the standard system of four equations || for 3 He nucleus in both the spin (S = 1/2) and 
isospin (T = 1/2) doublet states, as well as a separate system of two equations for the three-nucleon 
system with S = 1/2 and T = 3/2. 

The system of six spatial equations (fj) is the basic one for studying the 3 He nucleus (as well as the 
scattering and reactions in the system (2p, n) ) within the isospin formalism. In practice, however, 
precise calculations with the use of @ have some difficulties because of a considerable number of 
equations. We reformulate the problem in an equivalent form using the fact that the reason for an 
increase in the number of equations in the few-nucleon problem is the additional antisymmetrization 
of wave functions (the generalized Pauli principle for nucleons) connected with the introduction of 
isospin for treating a proton and a neutron as an isotopic doublet. It appears that not all the spatial 
components of the wave function in (fl|) are independent, and Eqs. (0) contain an implicit relation 
between them. The analysis of system (0) shows that it is suitable to introduce six new components 
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Then, as a result of the unitary (orthogonal) transformation, we get a new form of system @ with 
only two explicitly independent equations from the set of six ones. Respectively, only two spatial 
components of the wave function are independent, the rest being expressible through these two ones 
with the account of the permutations of coordinates: 

0'i (123) = (321) + ^V 3 ' (321) , ft (123) = ^V 3 (321) + (321) , 

02 (123) = -^3 (132) - ^V 3 ' (132) , $ (123) = -^0 3 (132) + i$ (132) . (5) 

We note an important fact that the unitary transformation (0) is of universal form and does not 
depend on the specific form of interaction potentials but only on the spin-isospin symmetry in ([[]), 
(0) . To find all the wave function components and to obtain the total antisymmetric function ([!]) for 
the system (2p,n), it is sufficient to solve the system of only two equations (instead of six ones) for 
two independent components 3 and 3 (further we introduce notations $1 = 3 and $2 = 0s)- 

On the other hand, the same functions $1 and $2 can be found directly from the system of two 
equations formulated in an equivalent representation which does not use the notion of isospin and 
treats a proton and a neutron as different particles. In this case, the total antisymmetric (in protons) 
function of the system (2p,n) with total spin 5 = 1/2 looks as follows: 



* (pi,P2,ra 3 ) = C'$i (pi,P2,n 3 ) + C"$2 (pi,P2,n 3 ) 



(6) 



where (' and (" are the two known components of the spin function of three nucleons with total 
spin S = 1/2. The spatial function <£>! (pi,P2,^3) is symmetric, while $2 (Pi,P2,^3) is antisymmetric 
with respect to permutations of identical protons. Within the representation free from the notion of 
isospin, the system of equations for $x and $2 has the following form: 
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In the system of equations (0) for 3 He nucleus, we take into account the dependence of charge- 
dependent nuclear potentials on spin and parity in the orbital momentum of nucleon pairs. In the 
case of 3 H nucleus consisting of two neutrons and one proton with the total spin 5" = 1/2, one has 
also Eqs. (0), but with exchanged indices n < — > p and without the Coulomb potential. Thus, Eqs. 
(0) completely determine the properties of three nucleons in the doublet spin state (S = 1/2) both 
in the bound states and scattering processes. 

The proposed approach is also applicable to the system of three nucleons in the quartet state 
(S = 3/2). All the properties of the system are determined in this case by only the triplet charge- 
dependent nuclear interaction potentials, and it is necessary to solve only one equation to find the 
spatial component of the wave function (the spin function of three nucleons is one-component in this 
case, and the spatial function of the system (2p, n) is antisymmetric with respect to permutations of 
the two protons). The equation looks as follows: 

$(123) + ^ J2 E^fn P )M [1±P(V)]*(123) = 0. (8) 

Z ij= 13,23+,- 

It should be noticed that the complete equivalence of the systems with different numbers of equations 
in different representations does not seem to be a miracle if one rewrites the nuclear interaction oper- 
ator in an equivalent form in terms of the Majorana exchange operator without explicit introduction 
of isospin operators. It follows from the general form of the central exchange interaction potentials 
presented in terms of the Majorana exchange operators that the transition to the equivalent repre- 
sentation without use of the isospin formalism is also possible for heavier nuclei and is reasonable for 
simplifying numerical calculations without any approximations. 

3 Results of calculations 

To study the properties of three-nucleon systems using the system of Eqs. (0), we develop an 
optimal variational scheme || with the use of a Gaussian basis. Precise calculations are carried out 
for the basic characteristics of the bound states of 3 H and 3 He nuclei for various nuclear interaction 
potentials. Table 1 contains the results of calculation of the binding energies B = —E and r.m.s. 
radii for the Afnan-Tang (ATS3) and Minnesota potentials (h 2 /M = 41.47 MeV- fm 2 in calculations 
with equal masses of a proton and a neutron, the Coulomb parameter being e 2 = 1.44 MeV- fm, 
and the parameters of potentials are taken from ||). The more realistic Afnan-Tang potential with 
interaction only in even orbital states is denoted by AT-(S3) + . The results for energies and radii are 
given with the accuracy of one unit in the last digit with high probability (the exact results for the 
binding energies are slightly greater, but within not more than 1 KeV). 

High accuracy of the calculations and the convergence of the results with spreading the variational 
basis are obtained with the use of a comparatively small number of basis functions (about 60-100 
ones for both the symmetric <3>iand antisymmetric $2 functions). It is found that the optimal way 
to achieve a given accuracy is to take into account about three times more Gaussian components for 
$1 than those for <3> 2 - Due to the approach based on system (|7|) without explicit use of isospin, our 
numerical results are of higher precision than the known Varga-Suzuki ones , although we used a 
less number of Gaussian variational basis functions. 

The greater the difference between the triplet and singlet interaction potentials and the greater 
the short-range repulsion, the more essential is the difference between the energies calculated exactly 
and those obtained in the spinless approximation. The complete calculation of energies with account 
of the difference in the neutron and proton masses (M p 7^ M n ) shows the binding energy of 3 H to 
increase a little and that of 3 He to slightly decrease, while r.m.s. radii are practically unchangeable. 
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For all used interaction potentials, the proton, neutron, and mass r.m.s. radii indicate the distinct 
specific structures of 3 H and 3 He nuclei. Namely, the proton radius in the 3 H nucleus is essentially 
smaller than the neutron one (something like the neutron "halo"), while the 3 He nucleus reveals the 
inverse pattern (the proton "halo"). The explanation of this effect lies in the fact that the attraction 
of a pair of nucleons in the singlet state is weaker than that in the triplet one. 

The obtained wave functions, having a suitable Gaussian representation, enabled us to calculate 
the main structure functions of the three-nucleon nuclei for various nuclear potentials. Fig. 1 presents 
the proton density distribution p p (r) as well as the neutron one p n (r) for a 3 H nucleus in the case 
of the AT-(S3) + potential (the nucleons are considered to be point-like particles). The peripheral 
neutron " halo" effect is seen distinctly due to the neutron density distribution being somewhat more 
long-range than the proton one (see also the r.m.s. radii in Table 1). At the same time, the proton 
density is essentially higher than the neutron one at the center of the nucleus. We have almost 
the same density distributions for the 3 He nucleus (even from the quantitative point of view), but 
with the exchanged protons and neutrons p <-> n. The main reason for such a regularity lies in the 
difference between the interactions in triplet and singlet states, while the Coulomb repulsion plays 
a negligible role. Fig. 2 shows the charge formfactors for 3 H and 3 He nuclei calculated with the 
same potential. Note that the formfactor profile of 3 He falls down already at small q 2 more rapidly 
as compared to that of 3 H, due to the greater radius of the proton density distribution in 3 He as 
compared with that in 3 H. And the formfactor of 3 He changes its sign at smaller q 2 (experimental 
value q^m, exp = 11.6 frrT 2 for 3 He) because of the essential role of short-range correlations between 
the protons. Fig. 3 presents the pair correlation functions of the 3 H nucleus both for the pair of 
neutrons ( nnij')) and for the neutron-proton pair ( g 2j np ij')) in the case of the AT-(S3) + potential. 
A noticeable difference between them is connected with a difference in the interactions in triplet and 
singlet states, and with the fact that, in the 3 H nucleus, the proton and a neutron interact mainly 
in the triplet state, while the two neutrons do in the singlet one. For comparison, we give also the 
deuteron wave function squared in the same figure for the same interaction potential (i.e., the pair 
correlation function of two nucleons in the triplet state; see the curve (np) D in Fig. 3) . The high 
accuracy of the known approximation based on two-particle correlation functions for the few-nucleon 
nuclei is confirmed once again. 

All the main conclusions are also valid for the Minnesota potential, in particular, those about the 
above-mentioned likeness of pair correlation functions of three- and two-particle systems. Moreover, 
for the both considered potentials, the pair correlation functions of 3 H are similar to those of 3 He 
(with the account of the substitution n <-> p like for density distributions). The structure of the nuclei 
is determined mainly by the triplet and singlet interaction potentials, while the role of the Coulomb 
interaction is small and leads to a little and almost proportional expansion of the 3 He nucleus in 
comparison with 3 H. But, in the case of the Minnesota potential having essentially less short-range 
repulsion than the AT-(S3) + one, the dips in the formfactors of 3 H and 3 He nuclei occur at greater 
q 2 , and the pair correlation functions have an essentially less decrease at short distances. 

4 Conclusions 

Thus, using the example of three-nucleon systems with central interaction potentials, we show the 
complete equivalence of the isospin formalism (with the total wave function being antisymmetric in 
the space of spin, isospin, and coordinate variables) and the proposed representation without explicit 
use of isospin (with the total wave function being antisymmetric in identical nucleons in spin and 
spatial coordinate variables). The obtained equations in the representation without use of isospin 
are unitary equivalent to the standard approach using the isospin representation. 
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The proposed approach is much more suitable (due to the small number of spatial equations 
that is determined by the dimension of a spin Young scheme) and enables one to carry out a precise 
study (with given accuracy) of three-nucleon systems in the bound states or scattering processes with 
various central nuclear interaction potentials. Using certain optimization schemes for the variational 
method with Gaussian basis, we carried out the precise calculations of the bound states of three 
nucleons with several nucleon-nucleon interaction potentials. The results are of higher precision than 
those available in the literature. Density distributions, formfactors, and pair correlation functions 
are calculated as well. 

The advantages of the proposed approach can reveal themselves in precise calculations of systems 
of four nucleons (when the total spin S = 0, there are two equations for spatial components), five 
nucleons (there are five equations at S — 1/2), and, perhaps, of six-nucleon systems (five equations 
for spatial components in the both cases of 6 He and 6 Li nuclei). The developed approach opens a 
real possibility to construct the realistic potentials describing completely, with reasonable accuracy, 
the main low-energy parameters of few-nucleon systems. 
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Figure captions 

Fig.l 

Profiles of proton and neutron density distributions for 3 H nucleus (for the AT-(S3) + interaction 
potential). Curves (1), (2) depict p(r), while curves (3), (4) show r 2 p(r). 
Fig.2 

Charge formfactors for 3 H and 3 He nuclei (for the AT-(S3) + interaction potential). 
Fig.3 

Pair correlation functions g2,np(f) and <?2,nn(?") for 3 H nucleus (for the AT-(S3) + interaction poten- 
tial). For comparison, the pair correlation function g2, PP (r) for 3 He and the deuteron wave function 
squared are shown. 
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Table 1 

Calculated binding energies B = —Eq (MeV) and r.m.s. radii (fm) for 3 H and 3 He nuclei. R p , 
R n , R m are the proton, neutron, and mass density distribution radii, respectively. (Experimental 
values: E = -8.481 MeV and R p = 1.57 fm for 3 H; E = -7.716 MeV and R p = 1.70 fm for 
3 He). 
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